On the violation of a local form of the Lieb-Oxford bound 
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In the framework of density-functional theory, several popular density functionals for exchange 
and correlation have been constructed to satisfy a local form of the Lieb-Oxford bound. In its 
original global expression, the bound represents a rigorous lower limit for the indirect Coulomb 
interaction energy. Here we employ exact-exchange calculations for the G2 test set to show that the 
local form of the bound is violated in an extensive range of both the dimensionless gradient and the 
average electron density. Hence, the results demonstrate the severity in the usage of the local form 
of the bound in functional development. On the other hand, our results suggest alternative ways to 
construct accurate density functionals for the exchange energy. 

PACS numbers: 71. 15. Mb, 31.15.cg 
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I. INTRODUCTION 

Density-functional theory [l[ (DFT) is one of the most 
popular electronic structure methods that, for a large 
variety of systems, produces accurate results with a rel- 
atively small computational cost. Functionals used for 
the exchange and correlation energy E xc [n] and potential 
V xc [n] play a central role in DFT 0. In fact, DFT be- 
came a mainstream method in Quantum Chemistry only 
after significant developments beyond the local-density 
approximation (LDA), first and foremost the generalized- 
gradient approximations (GGAs). 

Most GGAs are built in a way that -E^n] satisfies a 
set of known exact conditions. These conditions, together 
with an ansatz of the gradient form, enable GGAs to ac- 
count for inhomogeneities in the electronic density and 
thus improve upon the LDA in most cases. Two of the 
most prominent examples of the GGA are that of Becke 
from 1988 (B88) and the Perdew-Burke-Ernzerhof Q 
(PBE) functional. In the latter, seven known exact con- 
ditions to the S^cfn] are imposed to a simple ansatz. 
The result is an outstandingly accurate functional that 
performs well on a wide range of systems. Currently, 
PBE is the most popular functional for material appli- 
cations, whereas B88 - embedded in the B3LYP hybrid 
functional || - is the most popular functional in Quan- 
tum Chemistry 0. 

An important exact condition for S xc [n] is the Lieb- 
Oxford (LO) bound Q. The bound sets a rigorous 
lower limit for the indirect (quantum mechanical) part 
of the total Coulomb interaction energy. Hence, on one 
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hand, the LO bound is a fundamental condition in many- 
particle physics and relates to the analysis of the stability 
of matter [8j. On the other hand, the LO bound must 
be satisfied by all density functionals. In DFT this re- 
quirement can be conveniently formulated through the 
expression of the bound in terms of the LDA exchange 
energy (see below). 

Levy and Pcrdew 0, H| have suggested that in order 
to satisfy the LO bound, the exchange- energy density 
(exchange-energy integrand) must also be bounded. This 
local condition, however, is stronger than the global LO 
bound. In spite of the ambiguity, several functionals such 
as PBE rely on this local form to impose an extra con- 
straint on the functional. In fact, to the best of our 
knowledge, the LO bound has been applied solely in its 
local form in the development of functionals 0]. 

The use of the local LO bound has been questioned in 
several works [1(1 [H| 1 but the severity of the approxima- 
tion, i.e., the extent to which the local bound is violated, 
has not been explored until now. Addressing the simi- 
larities and differences between the global and local LO 
bound would be of particular importance in view of re- 
cent studies on the tightening of the LO bound [l2l fl3|. 
and the consequent developments or revisions of density 
functionals. 

In this work we use the well-known G2 test set [l4[ to 
analyze to which extent, and in which range of parame- 
ters, the local form of the LO bound is violated. The G2 
set of 148 molecules is commonly used to gauge the ac- 
curacy and predictive abilities of a given computational 
method, and it represents a broad range of chemical en- 
vironments [m, [H| . 

Our analysis is carried out in the following order. In 
Sec. [IT] we review and address the differences between the 
global and local forms of the LO bound within DFT. In 
Sec. lIIII we present our results for closed-shell molecules in 
the G2 set; here we solve exact exchange (EXX) energies 
within the Krieger-Li-Iafrate (l6| (KLI) approximation 
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and demonstrate the violation of the local bound. In the 
same context, we assess the local performance of differ- 
ent density functionals with respect to the enhancement 
factor. In addition, we compare the relation between the 
violation of the local bound and the distance from the 
global LO bound, and consider the spatial dependence 
of the violation of the local bound. Finally, in Sec. IIVI 
we summarize our findings and discuss how our results 
could be used in the development of density functionals. 



II. LIEB-OXFORD BOUND 

A. Global bound 

The LO bound in its original, global form applies 
to all three-dimensional (3D) non-relativistic, Coulomb- 
interacting systems. The bound can be expressed in 
terms of the indirect part of the interaction energy as 



or, alternatively, 



W XC [V] = h\V ee m - E H [n] >-C 



d 3 i 



,4/3 



(r) 



(1) 

where V ee = J2i>j l r « — r j\ 1 i s tne Coulombic interac- 
tion operator and Eh[it] is the classical Hartree energy. 
The expectation value of V ee is calculated over any nor- 
malized many-body wave function ^(ri, r/v) with the 
corresponding density n(r). For the prefactor C Lieb [l?} 
originally found a value C L = 8.52, but this was later re- 
fined by Lieb and Oxford H to C LO = 1.68 and numeri- 
cally by Chan and Handy Hi] to C* CH = 1.64. Recently, 
using nonrigorous but physical arguments the bound was 
tightened further to C = 1.44 (Ref. [3). 

Remarkably, the right-hand side of Eq. (p} has a form 
similar to the LDA for the the exchange energy, i.e., 



A / d 3 rn 4/3 (r) 



(2) 



with A = 3 4/3 7r _1 / 3 /4. Moreover, the left-hand side of 
Eq. (JTJl can be expressed as a density functional W^cfn] 

* ), so 



corresponding to the minimization of {^Sf T + V et 

that the ground-state density n(r) is produced. Following 
the definition of the exchange-correlation energy E xc we 
can now write 

W xc [n] < W xc [n] + T c [n] = E xc [n] < E x [n] , (3) 

where the first inequality is justified by the fact that the 
kinetic-energy part of the correlation energy is always 
non-negative, i.e., T c > 0. The second inequality fol- 
lows from the non-positiveness of the correlation energy 
E c [n] = E xc [n] — E x [n] < 0. This is straightforward to 
see in the constrained-search definition of the correlation 
functional fl9j |. 

Combining the above relations leads to a simple ex- 
pression of the global bound: 



E r 



E X DA [t 



< A, 



(5) 



where A = C/A. We note that, in this definition, A is 
a number that in principle is universal. Recent stud- 
ies [ill, [2(| on the LO bound have focused on finding 
maximum values for a density-functional A[n] or func- 
tion A(iV), with an aim to tighten the universal value for 
A. Indeed, the procedure in Ref. [l]| led to a conclusion 
that the bound can be tightened to A as 1.96 correspond- 
ing to C = 1.44. In the following, however, we will refer 
to the original LO value of C LO = 1.68 corresponding to 
A = 2.27. 



B. Local bound 

Before introducing the local bound it is useful to write 
the exchange-correlation energy in the form of the stan- 
dard GGA ansatz, i.e., 



E r 



d 6 rn{v)e x u c ^[n]F xc [n,Vn, 



(6) 



where e 3 {? EG is the exchange-correlation energy per par- 
ticle in the 3D homogeneous electron gas (3DEG), and 
JF XC is the enhancement factor including the corrections 
to the LDA. By definition it is non-negative, as well 
as its components T x and T c . In the standard GGA 
ansatz, T xc is written as a functional of n and Vn, or 
more conveniently, in terms of the dimcnsionless density 

1/3 

gradient s = \\7n\ / (2kpn) (with kp = (37r 2 n) be- 
ing the Fermi momentum) and the Wigner-Seitz radius 
r s = 3 1 / 3 (47rn) -1 / 3 . As indicated by the symbol ". . ." in 
Eq. ([B]), J- xc generally depends on other quantities and 
thus the expression refers to the exact E xc [n]. 

Combining the above integral form for E xc and the 
relations of the previous section we can rewrite the global 
LO bound as 



rfVn(r) e 3UJitj H Fx < A / tfVn(r) 



|e! DEG [nl 



(7) 

where we take the absolute values to deal with only non- 
negative quantities. The local bound suggested by Levy 
and Perdew 0, [H| (originally expressed for T xc ) is writ- 
ten for the integrands of Eq. ([7]), i.e., 



T x <\. 



(8) 



E x [n]>XE^ A [n] 



(4) 



If this local condition is satisfied, the global inequality in 
Eq. Q is trivially satisfied as well. However, the reverse 
implication obviously does not hold. In other words, 
Eq. © may be violated without the violation of the 
global bound. In particular, the local bound in Eq. dSJ 
might be considerably too strict to be a well-reasoned 
condition in the development of density functionals. This 
will be explicitly demonstrated in the following section. 
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FIG. 1: (Color online) Count histogram of exact enhance- 
ment factors for the electronic exchange as a function of the 
dimensionless gradient in the G2 test set. The results are 
compared to the enhancement factors of several density func- 
tionals. The line shows the local Lieb-Oxford bound that is 
violated for T x > 2.27. 



III. RESULTS 

We use the OCTOPUS code [2l[ to calculate the EXX 
energies and the corresponding EXX energy densities 
(per particle) e x xx for closed-shell molecules in the 
G2 test set. We use the norm-conserving Hartwigsen- 
Goedecker-Hutter pseudopotentials The EXX 

results are obtained from the optimized-effective- 
potential [23[ (OEP) scheme within the KLI approxi- 
mation that neglects the so-called orbital shifts in the 
full OEP. Apart from special cases such as long atomic 
chains [24[ the KLI approximation has been shown to be 
extremely accurate with respect to the full OEP [2|| . 

The obtained EXX energy density per particle is 
directly related to the exchange enhancement factor 



through e£ 



EXX 



C 3DEG 



T x . Hence, after calculating both 



the EXX and LDA results we can visualize the distribu- 
tion of the enhancement factors with respect to the local 
bound given in Eq. ((5]). The procedure is the following. 
For every coordinate r^, within every molecule, we col- 
lect the values r s (ri), s(rj), and J- X (^i), so that we can 
make a statistical count histogram of J- x as a function of 
s or r s . Figure [1] visualizes the situation as a function 
of s. The local LO bound is shown as a horizontal red 
(dark gray) line, so that all the values above the line at 
T x > 2.27 violate the local bound according to Eq. 
We find significant violation at s > 3. This range of the 
dimensionless gradient, estimating the ratio of the den- 
sity variation in the Fermi wavelength scale, corresponds 
to the tail of the electronic density. The local dependence 
is studied in more detail below. 

In Fig. Q] we also assess the enhancement factors of 
several popular density functionals with respect to the 
EXX results and to the local bound. The function- 
als include PBE (exchange only), B88, Gill's functional 
from 1996 [26| (G96), Wu's and Cohen's functional from 



FIG. 2: (Color online) Same as Fig.[T]but as a function of the 
Wigner-Seitz radius r s = 3 1/3 (47rn)~ 1/3 . 



2006 [27| (WC), and a Bayesian fit for the enhancement 
factor by Mortensen and co-workers [28| (BAY) . In Fig. Q] 
we can find a large variance in the behavior of different 
functionals, especially at large s that corresponds to the 
asymptotic exponential tail of the atomic or molecular 
charge distribution. 

As expected, PBE and WC obey the local bound by 
definition. However, they are relatively far from the his- 
togram of exact results when s > 2 and thus miss the 
correct asymptotic behavior of T x . BAY has a similar 
trend but follows the EXX to larger values of s, as it is a 
Pade fit to experimental atomization energies of a subset 
of the G2 test set. 

Interestingly, B88 performs best of all the tested ap- 
proximations. This is due to the fact that B88 has a pa- 
rameter chosen to reproduce the Hartrcc-Fock exchange 
energies for atomic systems that are conceptually similar 
to the G2 test set considered here. B88 also by con- 
struction obeys the correct asymptotic behavior of the 
exchange energy per particle, i.e., e x (r) — > — l/(2r) as 
r — > oo. G96 does not obey this limit and strongly devi- 
ates from the EXX results at large s. 

Figure [2] shows the distribution of the EXX enhance- 
ment factors as a function of the Wigner-Seitz radius r s . 
The local bound is violated at r s > 5 corresponding to 
densities around n < 0.002 bohr -1 . Considering a local 
version of a recently proposed tighter global bound [l3| 
with A = 1.96 (instead of the LO one) would lead to vi- 
olation already at r s ~ 4 (n ~ 0.004 bohr -1 ). In any 
case, the violation clearly seems to occur in the tail of 
the electronic density as shown explicitly below. 

In Fig. [3] we have a closer look at the enhancement fac- 
tor as a function of s for four diatomic molecules includ- 
ing H2, CI2, Na2, and Li2. H2 and CI2 violate the local 
LO bound more rapidly than Na2, and Li2. It is interest- 
ing to compare these results with the distance from the 
global bound considered by Odashima and Capelle [29l |. 
In their study H2 was found to be closest to the global 
bound with A[n] ~ 1.25. Therefore the high tendency of 
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FIG. 3: (Color online) Enhancement factor as a function of 
the dimensionless gradient s for four diatomic molecules. 
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FIG. 4: (Color online) Spatial behavior of the enhancement 
factor across the molecular axis of CO2. The density n and 
dimensionless gradient s are also shown (in arbitrary units). 
The arrows mark the position of the atoms. 

H 2 to violate the local bound in Fig.[3]is plausible and un- 
likely to be purely accidental, although obviously there is 
no rigorous implication from the global to local bound as 
discussed in Sec. [TTJ In fact, in Ref. [29[ L12 and CI2 were 
found to have A[n] « 1.21 and 1.11, respectively, whereas 
the local bound-breaking tendency is much higher in the 
latter system (Fig. [3]). Therefore, the spatial distribution 
of T x has an important role, so that there is no clear cor- 
relation in the sensitivity of a particular molecule to the 
global and local bounds. As discussed above, the only 
rigorous implication is the fact that if the local bound is 
obeyed (as in PBE), the global bound is obeyed as well. 

Finally we give an example of the spatial dependence 
of the enhancement factor. Figure [4] shows T x along the 



molecular axis of CO2, plotted together with the density 
and the dimensionless gradient (in arbitrary units) . The 
arrows mark the position of the atoms. Furthermore, for 
clarity we removed the points inside the cores of the pseu- 
dopotentials. The figure confirms our statement within 
Fig. [TJ that the local bound is violated in the tail of the 
electronic density where n is small and s is large. Sim- 
ilar spatial behavior was found in all the cases that we 
checked in this detail. 



IV. CONCLUSIONS AND OUTLOOK 

We have explicitly studied the conceptual difference 
between the universal Licb-Oxford bound and its local in- 
terpretation used in the development of density function- 
als. Our molecular examples in the G2 test set demon- 
strate that the local bound is broadly violated when exact 
exchange enhancement factors are considered. Therefore 
the use of the local bound in the development of func- 
tionals , e.g., in generalized-gradient approximations is 
questionable, even if the condition is straightforward to 
implement amd the fulfillment of the global Lieb-Oxford 
bound is guaranteed. In short, the local bound is simply 
much stricter that has been previously thought. 

Our exact-exchange results pinpoint the violation of 
the local bound to s > 3 or r s > 5 corresponding to the 
tail of the electronic density. We find no clear correlation 
between the distance from the global bound and the de- 
gree of violation of the local bound; this demonstrates the 
complexity of the enhancement factor in the tail region 
throughout the ensemble of molecules. 

The surprisingly uniform distribution of the exact en- 
hancement factors in Figs. [TJ and [2] suggest a construction 
of a density functional according to the observed J- X (s). 
This could be done with a fitting procedure or, better, 
by tailoring a physically motivated ansatz in the (meta-) 
generalized-gradient fashion that is able to reproduce the 
exact .Fx(s) to a reasonable extent. Such an ability is of 
particular importance when describing physical proper- 
ties that depend on the correct description of the tail of 
the electronic density, e.g., ionization and Rydberg exci- 
tations etc. Within the construction, however, a simul- 
taneous implementation of the global Licb-Oxford bound 
is a tedious (if not practically impossible) task. These 
aspects will be studied in detail in our future works. 
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